Decoherence of quantum registers 
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The dynamical evolution of a quantum register of arbitrary length coupled to an environment of 
arbitrary coherence length is predicted within a relevant model of decoherence. The results are 
reported for quantum bits (qubits) coupling individually to different environments ('independent 
decoherence') and qubits interacting collectively with the same reservoir ('collective decoherence'). 
In both cases, explicit decoherence functions are derived for any number of qubits. The decay of 
the coherences of the register is shown to strongly depend on the input states: we show that this 
sensitivity is a characteristic of both types of coupling (collective and independent) and not only of 
the collective coupling, as has been reported previously. A non-trivial behaviour ( "recoherence" ) is 
found in the decay of the off-diagonal elements of the reduced density matrix in the specific situation 
of independent decoherence. Our results lead to the identification of decoherence-free states in the 
collective decoherence limit. These states belong to subspaces of the system's Hilbert space that 
do not get entangled with the environment, making them ideal elements for the engineering of 
"noiseless" quantum codes. We also discuss the relations between decoherence of the quantum 
register and computational complexity based on the new dynamical results obtained for the register 
density matrix. 

PACS numbers: 03.65.Yz, 42.50.Lc, 03.67.-a, 03.67.Lx 



I. INTRODUCTION 

When any open quantum system, for example a quan- 
tum computer, interacts with an arbitrary surrounding 
environment, there are two main effects that have to 
be considered when examining the temporal evolution. 
First, there is an expected loss of energy of the initial 
system due to its relaxation, which happens at the rate 
T rei wnere T rei is the relaxation time-scale of the system. 
Second, there is a process that spoils the unitarity of the 
evolution, the so-called decoherence 0, where the con- 
tinuous interaction between the quantum computer and 
the environment leads to unwanted correlations between 
them in such a way that the computer loses its ability 
to interfere, giving rise to wrong outcomes when execut- 
ing conditional quantum dynamics. This phenomenon is 
characterized by a time that defines the loss of unitarity 
(i.e. the departure of coherence from unity), the deco- 
herence time Tdec- Usually, the time-scale for this effect 
to take place is much smaller than the one for relaxation, 
hence quantum computers are more sensitive to decoher- 
ence processes than to relaxation ones. For practical ap- 
plications in quantum computing there are several differ- 
ent systems that might provide a long enough r re i ; how- 
ever, what really matters for useful quantum information 
processing tasks (e.g. quantum algorithms) to be per- 
formed reliably is the ratio T gating /T dec . T gatingi the time 
required to execute an elementary quantum logic gate, 
must be much smaller than Tdec- As a rough estimate, 
fault-tolerant quantum computation has been shown to 
be possible if the ratio T gat i ng to Td ec of a single qubit is of 
the order of 10 -4 || (a qubit is a two-state quantum sys- 
tem, the basic memory cell of any quantum information 



processor). 

Decoherence and quantum theory are unavoidably con- 
nected. Indeed, the ubiquitous decoherence phenomenon 
has been ultimately associated with the "frontiers" be- 
tween the quantum behaviour of microscopic systems 
and the emergence of the classical behaviour observed in 
macroscopic objects roughly speaking, the decoher- 
ence time Tdec determines the energy and length scales 
at which quantum behaviour is observed. It generally 
depends non-trivially on several different factors such as 
temperature, dimensionality, quantum vacuum fluctua- 
tions, disorder, and others whose origin is less well known 
(hardware characteristics). The time-scale for decoher- 
ence depends on the kind of coupling between the system 
under consideration and the environment, in a range that 
can go from pico—seconds in excitonic systems || up to 
minutes in nuclear spin systems 

The discovery of algorithms for which a computer 
based on the principles of quantum mechanics [^) can 
beat any traditional computer, has triggered intense 
research into realistic controllable quantum systems. 
Among the main areas involved in this active research 
field are ion traps [|| , nuclear magnetic resonance (NMR) 
0, quantum electrodynamics cavities ||, Josephson 
junctions ||, and semiconductor quantum dots The 
main challenge that we face is to identify a physical sys- 
tem with appropriate internal dynamics and correspond- 
ing external driving forces which enables one to selec- 
tively manipulate quantum superpositions and entangle- 
ments. For this to be done, the candidate system should 
have a sufficiently high level of isolation from the en- 
vironment: quantum information processing will be a 
reality when optimal control of quantum coherence in 
noisy environments can be achieved. The various com- 
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munities typically rely on different hardware method- 
ologies, and so it is important to clarify the underlying 
physics and limits for each type of physical realization of 
quantum information processing systems. As mentioned 
above, these limits are mainly imposed by the decoher- 
ence time of each particular system. However theoretical 
work has shown that, in addition to quantum error cor- 
recting codes , there are two additional ways that 
may potentially lead to decoherence-free or decoherence- 
controllcd quantum information processing: one of them 
is the so-called "error avoiding" approach where, for a 
given decoherence mechanism (e.g. collective decoher- 
ence discussed below), the existence of "decoherence- 
free" subspaces within the system's Hilbert space can be 
exploited in order to obtain a quantum dynamics where 
the system is effectively decoupled from the environment 
[jTjLtni . This approach requires the system-environment 
coupling to possess certain symmetry properties. The 
other approach is based on "noise suppression" schemes, 
where the effects of unwanted system-bath interactions 
are dynamically controlled using a sequence of 'tailored 
external pulses' |l(||l7|]. These strategies have motivated 
much theoretical and experimental work. In particu- 
lar, there are some recent experimental results regarding 
the engineering of decoherence-free quantum memories 

us 

We devote this paper to the study of decoherence of 
an arbitrary quantum register (QR) of L qubits. In ad- 
dition to providing a general theoretical framework for 
studying decoherence, we examine in detail the two lim- 
its where the qubits are assumed to couple (i) indepen- 
dently and (ii) collectively to an external (bosonic) reser- 
voir. The reservoir is modelled by a continuum of har- 
monic modes. In Section II we show that the decoher- 
ence process is very sensitive to the input states of the 
register and give explicit expressions for the coherence 
decay functions. We have found that these functions 
possess a novel behaviour which we label "recoherence" 
(or coherence revivals) in the specific case of indepen- 
dent decoherence. This behaviour depends strongly on 
the temperature and the strength of the qubit-bath cou- 
pling. By contrast, for a certain choice of the QR input, 
the calculation of the reduced density matrix leads to the 
identification of decoherence-free quantum states [ p^|Jl5| 
when the qubits are coupled "collectively" to the envi- 
ronment, i.e. when the distance between qubits is much 
smaller than the bath coherence length and hence the 
environment couples in a permutation-invariant way to 
all qubits. The calculations in this paper were motivated 
by the results of Ref. [ p0[|2"l| ] . The present paper shows 
that some subtle but important details of these earlier 
results are incomplete. Particularly, the calculation of 
the L— QR density matrix reported here leads both to 
new qualitative results, when analyzing the behaviour of 
coherence decay, and new quantitative results, when es- 
timating typical decoherence times: these novel results 



emerge for L > 1, as discussed in Section III. Conclud- 
ing remarks are given in Section IV. We emphasize that 
the results of this paper are not restricted to a particular 
physical system; they are valid for any choice of the qubit 
system (e.g. photons, atoms, nuclei with spin 1/2, etc.) 
and any bosonic reservoir. 



II. INDEPENDENT VERSUS COLLECTIVE 
DECOHERENCE 

Let us consider the general case of a L— QR coupled to 
a quantized environment modelled as a continuum of field 
modes with corresponding creation (annihilation) bath 
operator w (b). Throughout this work we will analyze 
pure dephasing mechanisms only. We will not consider 
relaxation mechanisms where the QR exchanges energy 
with the environment leading to bit-flip errors. Hence, 
the n th qubit operator J n = J™ (J™ = J™ = 0) @. 
As mentioned earlier, this is a reasonable approach since 
decoherence typically occurs on a much faster time scale 
than relaxation. The dynamics of the qubits and the 
environment can be described by a simplified version of 
the widely studied spin-boson Hamiltonian p3| : 



H = ^ 6 nJz 
n=l 



where the first two terms describe the free evolution of 
the qubits and the environment, and the third term ac- 
counts for the interaction between them. Here q n de- 
notes the coupling between the n th qubit and field modes, 
which in general depends on the physical system under 
consideration. The initial state of the whole system is 
assumed to be p s (0) = p®(0) ® p B (0) (the superscripts 
stand for system, qubits, and bath) , i.e. the QR and 
the bath are initially decoupled p4j| . We also assume 
that the environment is initially in thermal equilibrium, 
a condition that can be expressed as 



Ao)=n^( r )=n 



exp(-f3huj k bjb h 



(2) 



where (iV w ) = [ exp(— /3hu> k ) — l] 1 is the Bose-Einstein 
mean occupation number, ks is the Boltzmann constant 
and (3 = 1/fceT. For the model of decoherence presented 
here it is clear that we are in a situation where the qubit 
operator J™ commutes with the total Hamiltonian H and 
therefore there is no exchange of energy between qubits 
and environment, as expected. We will not attempt to 
perform a group-theoretic description for the study of 
quantum noise control (25). Instead, we study the specific 
real-time dynamics of the decay in the QR-density matrix 
coherences within the model given by the Hamiltonian 
Eq. (1) (i.e. Abelian noise processes, in the language of 
Refs. |6p7)). 
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In the interaction picture, the quantum state of the 
combined (qubits + bath) system at time t is given 
by = U I {t)YH{Q)) 1 where |*(0)) is the initial 

state of the system, L/ f (t) is the time evolution oper- 
ator, Uj(t) = Texp [ - i/h f Q H^t^dt'], and T is the 
time ordering operator. For the Hamiltonian (1) we 
introduce the notation H = Ho + V, where Hq = 
J2n=i e nJz + J2k € >eblb k denotes the free evolution term 

and V — J2 n k Jz(.9h"k 3fc*^ fc ) ^ s tne interaction term. 
Hence, the interaction picture Hamiltonian is given by 
H,{t) = ul{t)VU {t), with U a = &qp(-%H t). A simple 
calculation gives the result 



(3) 



which allows us to calculate the time evolution operator 
Uj(t). The result is (see appendix B.l for details) 



U I (t)=exp i$ u ft) 



exp 



E J "{e(*)&£-em} 



(4) 



with 



<M') = EI J X 



,2 t^t - sin(^) 
1 - e""** 



(5) 



This result differs from the one reported in Ref. pQ ] 
where the izme ordering operation for U t (t) was not per- 
formed. As will become clear later, this correction alters 
the resulting calculation of Ref. |^0| , and hence changes 
the results for the reduced density matrix of an arbitrary 
L— QR. Based on the time evolution operator of Eq. (4), 
we report here a different result for this density matrix 
and discuss its implications with respect to those of Ref. 

Unless we specify the contrary, we will assume that 
the coupling coefficients g^ (n = 1, ...,£) are position- 
dependent, i.e. that each qubit couples individually to 
a different heat bath, hence the QR decoheres 'indepen- 
dently' . Implications for the 'collective' decoherence case 
will be discussed later. Let us assume that the n th qubit 
experiences a coupling to the reservoir characterized by 



9'L 



h exp( 



-ik-r r 



(6) 



where r n denotes the position of the n th qubit. It is 
easy to see that the unitary evolution operator given by 
Eq. (4) produces entanglement between register states 
and environment states ^Oj. The degree of the entan- 
glement produced will strongly depend on the qubit in- 
put states and also on the separation ||T" m — r n \\ between 
qubits because of the position dependent coupling. As 



we will show later, for some kind of input states no de- 
coherence occurs at all despite the fact that all of the 
qubits are interacting with the environment. We will 
also identify states with dynamics decoupled from ther- 
mal fluctuations; this fact may be relevant when design- 
ing experiments where the involved quantum states have 
dephasing times (mainly due to temperature dependent 
effects) on a very short time scale, as in the solid state 
for example. We will also see that the above features 
are key issues when proposing schemes for maintaining 
coherence in quantum computers |H| . 

Due to the pure dephasing (i.e. Abelian) characteristic 
of the noise model considered here, we can calculate an- 
alytically the functional dependence of the decay of the 
coherences of the QR by taking into account all the field 
modes of the quantized environment. We shall follow the 
notation of Ref. p3] . Let us consider the reduced density 
matrix of the L— QR: the matrix elements of this reduced 
operator can be expressed as 



,h\Tr B {p S (t)}\j L ,J L 



,Ji) 



(7) 



where {i n ,jn} = {h,ji;i2,j2;—\i t ,j L } refers to the 
qubit states of the QR and 



p s (t) = U I (t)pQ(0)^p B (0)U}(t) 



(8) 



From this expression it becomes clear that the dynamics 
of the register is completely determined by the evolution 
operator U 2 (t). In Eq. (8), the initial density matrix 
of the register can be expressed as p® (0) = pf ■ (0) ® 

1 '-'1 

(0)®...<g>p? ,• (0), where pf ■ = \i n ){j n \. In this 
= | ± h) are the possible eigenstates of 
J™ and will be associated with the two level qubit states 
|1) and |0), respectively. The eigenvalues i n = ±-| and 
j n = ±i. In what follows, the subscripts of Eq. (7) will 
be renamed with the values 1 and to indicate the actual 
values I and —4, respectively. Hence we can rewrite Eq. 
(7) as 



Pi 2 <H 

expansion, \i 



(t) = Tr E 



p B (0)U^(t)Uf^(t) 



Q 



(0) 



(9) 



t/K > (t) = exp [i £ \g k \ 2 s(u; k , t) £ e^^i, 

fc m.n 

exp [J2 [g k ^ k (t)i n e- ik - r "bt - g yZ k {t)i n j k ^b h ) 



(10) 

and calculate explicitly the decay of the coherences for 
the L— QR. The result is (see appendix B.2) 



P? in , 3 „ } (*) = ex P - E M C K'*) 



(11) 
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coth 



/ Tiui k \ , , . . . . . , 
exp ^ ^ \g k \ 2 s(Lu k ,t)(i m i n - j m jn) cos k ■ r 

fc;m,n 

exp[-2i ^2 \g h \ 2 c(uj k ,t)i m j n smk ■ r mn P® nijn} {0), 



where r mn = r m —r n is the relative distance between the 
m th and n th qubits, s(uj k ,t) — [uj k t - sm(uj k t)]/ (huj k ) 2 , 
and c(uj k ,t) = [1 — cos(w fe i)]/(?i(jj fc ) 2 . In the continuum 
limit, Eq. (11) reads 



pg nJn} (i) = exp i{6 d (M a )- A d (i,i»)} 



exp 



-r d (t,t s ;T) 



(12) 



where 

<d d (t,t s ) = I dxF d {* •■ 



2 ^ ^ {imin jmjri) COs( 



(13) 



A d (t,t s ) = 2 dujl d {uj)c(u,t) ^ imjnsin(uit s ) , (14) 

n = 2 



and 



r d (t,t a ;T)= / dw/ d (w)c(a;,t)coth(^-)x (15) 



2 (i m - jm)(in ~ jn) COs(wt s )|. 



X=l,Tl=2 



Here u; T = k B T/fi is the thermal frequency (see discus- 
sion below), Lot s = k ■ r mn sets the "transit time" t s , 
and = E fc ^-^)l5j 2 = f GM|. 9 H| 2 is the 

spectral density of the bath. This function is charac- 
terized by a cut-off frequency uj c that depends on the 
particular physical system under consideration and sets 
I d {oj) i— » for oj >> u c H|. We see that an explicit 
calculation for the decay of the coherences requires the 
knowledge of the spectral density I d (to). Here we assume 
that I d {uj) = a d oj d e~ u) l u)c j22|, where d is the dimension- 
ality of the field and a d > is a proportionality constant 
that characterizes the strength of the system-bath cou- 
pling. Hence, the functional dependence of the spectral 
density relies on the dimensionality of the frequency de- 
pendence of the density of states G{lo) and of the coupling 
g(uj). In Eq. (12) the "transit time" t s has been intro- 
duced in order to express the QR-bath coupling in the 
frequency domain. This transit time is of particular im- 
portance when describing the specific "independent" de- 
coherence mechanism, because of the position-dependent 
coupling between qubits and bath modes. However, in 
a scenario where the qubits are coupled "collectivelly" 



to the same environment, t s does not play any role (see 
below). 

The result of Eq. (12) differs in several respects 
from the one reported in pc|| : the decoherence func- 
tion T d (t,t s ]T) contains additional information about 
the characteristics of the independent decoherence and 
the way in which the individual qubits couple to the en- 
vironment through the position-dependent terms which 
are proportional to cos(fc • r mn ). In essence this means 
that the entanglement of the register with the noise 
field depends on the qubit separation. Also the ex- 
pression for pQ . } (t) reveals the new dynamical factor 
tt d (t,t s ) = <d d (t,t s ) — A d (t 7 t s ) which must be taken into 
account when determining typical decoherence times for 
the L-QR. 

It is interesting that the decoherence effects aris- 
ing from thermal noise can be separated from the 
ones due purely to vacuum fluctuations. This 
is simply because the average number of field ex- 
citations at temperature T can be rewritten as 
{Nj) T = iexp(-ftw/2fc s T)cosech(7kj/2fc B T), and 
hence coth(?iw/2fc B T) = 1 + 2(N W ) T in Eq. (15). The 
other term contributing to the decay of the coherences 
in Eq. (12), Nd(t, i s ), is due purely to quantum vacuum 
fluctuations. The separation made above allows us to 
examine the different time-scales present in the (QR + 
bath) system's dynamics. The fastest time scale of the 
environment is determined by the cut-off: r c ~ u>~ , i.e. 
r c sets the "memory" time of the environment. Hence, 
the vacuum fluctuations will contribute to the dephasing 
process only for times t > t c . Also note that the charac- 
teristic thermal frequency lo t = ksT/h sets another fun- 



damental time scale r„ 



: thermal effects will affect 



the qubit dynamics only for t > t t . Hence we see that 
quantum vacuum fluctuations contribute to the dephas- 
ing process only in the regime r c < t < t t . From this 
identification it becomes clear that the qubit dynamics 
and hence the decoherence process of our open quantum 
system will depend on the ratio of the temperature-to- 
cut-off parameter lo t /uj c and the spectral function I d {u>). 
Later we will analyze how different qualitative behaviours 
are obtained for the decoherence depending on the rela- 
tionship between the cut-off and the thermal frequency. 
It is worth noticing that the analytical separation be- 
tween the thermal and vacuum contributions to the over- 
all decoherence process is a convenient property of the 
pure dephasing (Abelian) model considered here. Fur- 
ther generalizations of this model, e.g. by including re- 
laxation mechanisms, should make this separation non- 
trivial because the problem becomes no longer exactly 
solvable. 

Next, we analyze the case of "collective decoherence". 
This situation can be thought of as a bath of "long" co- 
herence length (mean effective wave length A) if com- 
pared with the separation r mn between qubits, in such a 
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way that A >> r mn and hence the product of Eq. (6) 
has exp(— ik ■ r n ) ~ 1. Roughly speaking, we are in a sit- 
uation where all the qubits "feel" the same environment, 
i.e. g? = g k . A similar calculation to the one followed 
in Appendix B gives the following result for the decay of 
the coherences for the case of collective coupling to the 
reservoir: 



pf (t) = exp <^ iQ d (t) 



E 



in 

771 — 1 

r L 



expi -T d (t;T) 



where 



E 1 

m— 1 



G d (t) = / dL0l d {ij)s{L0,t) , 



m= 
jm) 



and 



r d (t;T) = J <L>I d (uj)c(u),t)coth(^- 



(16) 



(17) 



(18) 



Expressions (12) and (16) are to be compared with those 
reported in Refs. [50] and Clearly, the new result 

for the evolution operator Ui(t) of Eq. (4) induces a 
QR-environment dynamics different to the one reported 
in ^ZQ. This fact has been pointed out in Ref. |2l[] for 
the situation of collective decoherence; in this particu- 
lar case, our results coincide with theirs. However, we 
obtain different results when considering the situation 
of independent decoherence: the new dynamical factor 
^d(t, t s ) includes additional information about individual 
qubit dynamics that cannot be neglected. Indeed, from 
the results derived in Ref. pl| , the authors argue that 
the damping of the independent decoherence is insensi- 
tive to the type of initial states and hence the sensitivity 
to the input states is only a property of the collective 
decoherence. As can be deduced from Eqs. (12) and 
(16), we find that this statement is not generally true 
and that the sensitivity to the initial states is a property 
of both collective and independent decoherence. This re- 
sult is particularly illustrated for the case of a 2— QR in 
the next section. From the expressions (17) and (18) we 
see that for hu> << ksT, the high-temperature environ- 
ment (high-TE), the phase damping factor T d (t; T) is the 
main agent responsible for the qubits' decoherence while 
the other dynamical damping factor Q d (t) plays a minor 
role. In this case oj c is actually the only characteristic 
frequency accessible to the system (u) c < < oj t ) and ther- 
mal fluctuations always dominate over the vacuum ones. 
However, when we consider the situation lj c >> lu t , the 
low-temperature environment (low-TE), these damping 
factors compete with each other over the same time scale, 
and £>d(t) now plays a major role in eroding the qubits' 
coherence. Here there is a much more interesting dynam- 
ics between thermal and vacuum contributions (see next 



section). This shows the difference and the importance 
of the additional terms of the reduced density matrix re- 
ported here when compared with those of Refs. po| , pT[ . 
The above statements will be illustrated in the next sec- 
tion. 

So far, the dynamics of the qubits and their coupling 
to the environment has been discussed in the interaction 
representation. To go back to the Schrodinger represen- 
tation, recall p Sch {t) — t/ (i)y0 7 (i)tTj'(i), with p I (t) as cal- 
culated before for the qubits decoherence (Eqs. (12) and 
(16)). Also note that in the Schrodinger representation 
(here denoted with subindex Sch) U (t) introduces mix- 
ing but not decoherence. Next, we consider some par- 
ticular cases which allow us to evaluate the expressions 
(12) and (16) and hence give a qualitative picture of the 
respective decoherence rates for both collective and inde- 
pendent decoherence situations. 



III. DIMENSIONALITY OF THE FIELD AND 
DECOHERENCE RATES FOR FEW QUBIT 
SYSTEMS 

In this section we analyze the qualitative behaviour of 
the decay of the coherences given by Eqs. (12) and (16) 
for single and two qubit systems. 



A. Single qubit case 

Here we consider the case of only one qubit in the 
presence of a thermal reservoir, as defined in Eq. (2). 
Hence, for both types of coupling (Eqs. (12) and (16)) 
we get 



plJt)=eM-r d (t,T)}pQjO) 



(19) 



By using Eq. (4) , with d = 1, it is easy to show that the 
populations remain unaffected: p®(t) — p®(0), i = 0, 1. 
In general (i ^ j), the decay of coherences is determined 
by Eq. (19). Here we can identify the main time regimes 
of decoherence for different dimensionalities of the field. 
In what follows, we consider the case of reservoirs with 
both one-dimensional density of states ("Ohmic") and 
three-dimensional density of states ("super-Ohmic"), i.e. 
G(uj) — constant and G(tu) = lu 2 , respectively, where the 
frequency dependent coupling g{ui) oc ^Juj, as considered 
in [^o). From Eq. (19) we see that a general solution 
for the case d = 1 requires numerical integration (see 
Appendix A). However, for the case where the interplay 
between thermal and vacuum effects is more complex, 
i.e. the low-TE (lo t << w c ), we can solve it analytically. 
Here we get 



ri(i, T)« Cl 



2uJ T t arctan(2oj T i) + i In 



1 + UJ 2 c t 2 

l + Au 2 t 2 



(20) 
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where ci = cti/Ti 2 . On the other hand, an exact solution 
for the super-Ohmic case d = 3 (Eq. (19)) can be found 
for any temperature T. The result is: 



r 3 (t, T)=cJ 2 C(2, 0) + C(2, 1 + 0) - C(2, 6 + iu T t) 



((2,9-iuj T t) 



1 - urit 



2+2 



H 2 ] 



(21) 



where C(«,g) = [1/T(z)] J™ dt [^-^"^/(l-e"*)] is the 
generalized Riemann zeta function, T(z) = J o °° cfa f z_1 e _t 
is the Gamma function, and C3 = a^u) 2 /h 2 . 

For the purpose of this paper, we concentrate on the 
case L > 1 for which we have several novel results. We 
leave the analysis of the L = 1 case to Appendix A, where 
we discuss the process of identifying typical decoherence 
times for a single qubit, and the interplay between the 
different decoherence regimes as a function of the tem- 
perature. 



B. L = 2 qubit register 



(hi) i a ^ j a , i b ^ j b 
where h 



*aja:hji, = 2[l + (i - j a )(ib - jb) cosk ■ r a b\- 
Clearly, hw,io = ^01,01 = 2[1 + cosfc • r ab ], and 
= hoi,w = 2[l-cosfc-r a6 ]. Hence pf ajaiibjb (*) 
also shows collective decay. 

In the above cases, analytic expressions for the corre- 
sponding decoherence functions can be found. As in Sec- 
tion III (a), we shall consider two different surrounding 
environments. In the low-TE regime, the 'Ohmic envi- 
ronment' (d = 1) induces the following coherence decay 
(the high-TE requires numerical integration): 



p± t (t) «exp 



Ti(t, T) ± ic\ ( — arctan(o; c t_) 



— arctan(ttu£ . 
2 K ^ 



UJrt 



1 + ^2 



PL(P) ( 22 ) 



for i a = j a , and if, 7^ j^. For brevity, we have dropped 
the subindices of the reduced density matrix, and set 
t, — t s + 1, and t_ = t s — t. For i a ^ j a , and ib ^ jb we 
obtain the result 



Let us analyze the case of two qubits in the presence 
of the bosonic reservoir discussed in the present paper. 
Using the same expressions for the density of states G(ui) 
and for the frequency dependent coupling g{uS) as above, 
we will analyze the coherence decay for several different 
input states. We set the qubits at positions r a and r b with 
coupling factors given by <?" = g k e~ l ^' rn , and n = a,b. It 
is easy to see that the unitary evolution operator induces 
entanglement between qubit states and field states: Ui(t) 
acts as a conditional displacement operator for the field 
with a displacement amplitude depending on both qubits 
of the QR, as discussed in more detail in . As we have 
pointed out previously, it is this entanglement that is re- 
sponsible for the decoherence processes described in the 
present paper. In particular, the case of two qubits has 
P? a ia,i b i b ^ = ( i a^h\ Tr B {p S { t )}\j a d b )- Next we analyze 
the register dynamics for the limiting decoherence situ- 
ations described above. First, let us study the case of 
independent decoherence: 

(i) ia = ja, ib ^ jb- 

P? aia ,n ib ® = exv[iQd(t)fi aia , ibjb - T d (t;T)] x 
P? aiaAbjb (Q)i where f iaia .i bjb = 2i a (i b -j b ) cosk-r ab . 
Hence / 00 ,oi = /n,io = cosfc • r ab , and / 00 ,io = 
/n,oi = — cosfe • r ab : pf . . . (t) shows collective 
decay. This result is contrary to the one reported 
in Ref. fl20tl. 



(ii) i a 
p Q 



Ja, lb -- 



jb- 

-. nQ 



(0); as expected, the popula- 



tions remain unaffected. 



p± t (t) «exp 



2T 1 (t,T) T 
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111 



1 + u 2 c il 



In 



1 + u\t\ 

1 + 4^2 



, ^2t s arctan(2cj T t s ) — t_ arctan(2w r i_) — 
rctan^ij}) p± =1 (0) . (23) 



t, ar 



In Figs. 1 and 2 we have plotted the decay of two 
qubit coherence due to the coupling to an environment 
of the Ohmic type (Eqs. (23)). Here, Tf (t, T) are defined 
fromEqs. (23) asp± =1 (t) = exp[-Tf (t,T)]p+ =1 (0). From 
these figures we can see the variations of the onset of de- 
cay when increasing the temperature. Figure (1) shows 
how the departure of coherence from unity changes with 
t s (plot 1 (i)) whilst for high temperature (plot 1 (hi)) 
there is no variation with t s at all. In the limit of large t s 
(see Table 1), we recover the onset of decay of Fig. 6 (Ap- 
pendix A). Note the difference between the time scales of 
plots 1 (i) and (iii), and how an estimation of typical de- 
coherence times strongly relies on the temperature of the 
environment. Figure (2) shows how the coherence decay 
shown in Fig. 1 disappears for small t s values, i.e. for a 
given temperature, it is possible to find a t s from which 
there is no decoherence of the two qubit system. This 
interesting behaviour occurs only for the density matrix 
elements (10|Tr B {p s (f)}|01) = (01|Tr B {p s (t)}|10). 

In Table 1 we give some typical two qubit decoherence 
times Tdec for an Ohmic environment in terms of the tcm- 



G 



perature, coupling constants ci, and t s . In all of the ta- 
bles in this paper, we have taken the beginning of the de- 
coherence process to occur when the reduced density ma- 
trix of the whole system exhibits a 2% deviation from the 
initial condition, i.e. when exp[— Ti(t = Td ec ,T)] = 0.98. 
The tables have been generated from the corresponding 
decoherence functions reported in this paper. Here, tf is 
defined from exp[-r\(i = t f ,T)} = 0.01, i.e. the differ- 
ence between t / and Tdec gives an estimate for the dura- 
tion of the decoherence process, say tdecay] an d T^ ec , and 
are evaluated from the respective decoherence func- 
tions Vf(t, T), with i = 1,3. In order to gain insight into 
some characteristic time-scales, consider for example the 
case of the solid state, where in many situations the noise 
field can be identified with the phonon field. Here the 
cut-off lu c can be immediately associated with the Debye 
frequency uo D . A typical Debye temperature 6_d = 100 K 



has 



10 



13 



so e 



J T /Uc 



10~ 2 T. Hence 



we can see from Table 1 (Figs. 1 and 2) that for c\ = 0.25, 
uo c t s = 0.5, and T = 0.1 K, the decoherence process 
starts at t£ ps 23.5 fs and t7 w 43.7 fs, and it lasts for 



t 



decay 



10.3 ps (for t decay , expt-I 1 ! (t 



,0.1 K)] sat- 



urates above 0.01). If the strength of the coupling goes 

+ c «0.2ps, and t+ ecay « 3.5 
and t~7„„„„. are not reported 



down to c\ = 0.01 then r 
ns. In this latter case 



dec' decay 

since the coherence saturates to a value above 0.98. From 
the data reported in Table 1 is clear that the weaker the 
coupling between the qubits and the environment, the 
longer the decoherence times and the slower the decoher- 
ence process. In addition the higher the temperature, the 
faster the qubits decohere, as expected intuitively 

The effect of the transit time becomes more evident 
from Table 1: for large t s values there is no difference 
between r+ c and r dec (also i+ w t d ). This is be- 



cause for large t s , the contribution due to terms involving 
t s in Eq. (23) is negligible, hence the reduced density ma- 
trix pf_ 1 (t) = p d=1 (t) ~ exp[— 2Ti(t, T)] and hence has 
a similar behaviour to the single qubit case (Eq. (19)). 
Hence for large t s (e.g. t s = 10 4 ), Figs. 1 and 2 resem- 
ble the onset of decay of the single qubit case (see later 
Fig. 6) where there is no dependence on t s . We note 
that exp[— r+(t, T] (exp[—TT(t,T]) is the corresponding 
decay of the coherences p lo w = p mm (p lo m = p 01A0 ), 
hence P? o w (t) = p? 0i01 (t) for large t s . 

By contrast, the super-Ohmic d = 3 field leads to the 
following coherence decay 



P d =s(*) = eX P 



IC 3 



-r 3 (t,r)± 

sin(2 arctan(w c i_)) sin(2 arctan(w c i + )) 



2(l+o# 2 _) 
2u) c t cos(3 arctan(w c t s )) 
(l + ^?) 3/2 
for i a = j a , and i b ^ j b , and 



2(1 + ^) 



(24) 
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0.5 


0.436919 


saturates 


0.235446 


103.507 


0.25 


10° 


0.5 


0.183755 


saturates 


0.104119 


2.05958 


0.25 


10~ 3 


10 4 


0.290113 


1279.63 


0.290113 


1279.64 


0.25 


10° 


10 4 


0.127778 


3.45901 


0.127778 


3.45901 


0.1 


10~ 3 


0.5 


0.913573 


saturates 


0.37654 


2025.75 


0.1 


10° 


0.5 


0.303135 


saturates 


0.16504 


4.28334 


0.1 


10~ 3 


10 4 


0.47316 


5669.66 


0.473159 


5670.15 


0.1 


10 


10 


0.203549 


T 0/"* f n/~> 

7.86596 


0.203549 


7.86596 


0.01 


10~ 3 


0.5 


saturates 


saturates 


1.45274 


35004.7 


0.01 


10° 


0.5 


saturates 


saturates 


0.538502 


37.2732 


0.01 


io- 3 


10 4 


2.55738 


saturates 


2.55738 


40816.8 


0.01 


10° 


10 4 


0.709492 


73.8325 


0.709492 


73.8325 



Table 1. Characteristic times for two-qubit independent decoher- 
ence T^, c for d = 1 dimensional density of states of the field. Dif- 
ferent temperature, transit time, and coupling strength values are 
considered. i a ^ j a , ib ^ jb (see text). 



pf =3 (t) = exp 



2r 3 (t,r)T2c 3 



l - wit 



2+2 



+ 



+ 



2[1+C^] 2 2[l+^ 2 _] 2 
Q2 ( 

— 1 2C(2, 9 - iu T t s ) + 2((2, 9 + iLu T t s ) - 

C(2, e + iuj T t + ) - c(2, e - ito T t + ) - 



<(2, 9 + iuJ T t-) - C(2, 9 - ioj T t-) 



Pts(0) 

(25) 



for i a 7^ j a , and i b ^ j b . The results of Eqs. (24) and (25) 
are exact: no approximations have been made in obtain- 
ing them. Therefore, they are valid for any temperature 
of the environment. 

In Figs. 3 and 4 we have plotted the decay of two 
qubit coherence due to the coupling to a reservoir with 
three-dimensional density of states (Eqs. (25)). Some 
estimates for r c and tdecay have been given in Table 
2. The decoherence functions Tf(t,T) are defined as 

P f =3 (t) = eM-rf(t,T)]pt 1 (o). " 

As we can see, there is a non-monotonic behaviour for 
the decay of coherence for low temperature values, as 
can be seen from Figs. 3 (i), and (ii). The decay given 
by the functions T+(t,T) (Fig. 3) and T^(t,T) (Fig. 
4) saturates to a particular value, which is fixed by the 
strength of the coupling and the temperature of the reser- 
voir: the lower the temperature, the slower the decay, 
and the higher the residual coherence. Some estimates 
for these saturation values e~ Fs are shown in Table 

2. From Fig. 4 it is possible to find small t s values for 
which the onset of decay does not change in time and co- 
herence remains unaffected. This result is very different 
from the case of Fig. 3, where coherence either vanishes 
(at high temperatures) or saturates to a residual coher- 
ence value (at low temperatures). Also note that whilst 
nothing happens to the onset of decay of Fig. 4, the 
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coherence decay is amplified in the case of Fig. 3, for 
small t s values. From Figs. 3 and 4 (Figs, (i), and 
(ii)) we see that there is saturation of the decay in the 
presence of a non-trivial coherence process. However, for 
high temperatures (Figs. 3 (iii) and 4 (iii)) there is a 
monotonic behaviour where no saturation occurs at all 
and the residual coherence vanishes. 
We note that the non-monotonic behaviour reported here 
for the low temperature regime is not just a characteristic 
of high dimensionality fields: it also occurs for the d = 1 
field, as can be seen from Figs. 1 (i) and 2 (i). We believe 
that this purely quantum-mechanical phenomenon is as- 
sociated with an interplay between the quantum vacuum 
and thermal fluctuations of the system and the quantum 
character of the field. The system undergoes a dynam- 
ics where the environment manages to 'hit back' at the 
qubits in such a way that the coherences then exhibit an 
effective revival before vanishing (at high temperatures) 
or saturating to a residual value (at low temperatures). 
An essential feature of the model studied here is the fact 
that the QR and bath are assumed to be initially uncor- 
rected. In future work we will analyze the behaviour of 
these 'recoherences' with more general initial conditions, 
where the initial state of the combined system is allowed 
to contain some correlations between the bath and the 
QR (see also Ref. |2^]). Such studies will help to clar- 
ify the origins of this 'recoherence effect' and also the 
effectiveness of decoherence as a function of the initial 
conditions. We note that there is a previous work by Hu 
et al. |30| ] where the study of quantum brownian motion 
in a general environment gives rise to a similar behaviour 
(although in a different context) to the one reported here 
for the 'recoherences'. A more detailed analysis of the 
physical implications of this interesting behaviour of the 
coherence decay is intended to be addressed elsewhere 

It is interesting that the dynamics of coherences de- 
pends so strongly on the strength of the coupling. It can 
be seen from Table 2 and Fig. 5 that i) the coherences 
saturate (sat.) to a very high 'residual value' and show 
less than a 2% decay for weak coupling (e.g. C3 = 0.01) 
and low temperatures (T = 0.1 K). This result is to be 
compared with the d — 1 case where the coherences al- 
ways vanish for long i s ; ii) even at high temperatures 
(T = 10 4 K) where the coherences vanish, and weak cou- 
pling, we find the appearence of the 'recoherence effect' 
discussed above. This makes more evident the role of the 
QR-bath coupling: recoherences should appear only un- 
der certain conditions imposed by the strength C3 (hence 
by the spectral density) and the temperature. These con- 
ditions can be directly obtained from the decoherence 
functions reported in this paper. Typical T^ ec times for 
this d — 3 dimensional environment are given in Table 2. 
From these we can conclude that the two QR decoher- 
ence time scales are shorter, and the decoherence process 
occurs faster, than in the single qubit case. 
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sat. 
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10 2 
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0.066994 
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0.07645 


sat. 


0.449 


0.01 


10~3 


10 2 


9.7767 


sat. 


0.9802 


9.7767 


sat. 


0.9802 


0.01 


10 2 


10 2 


0.07124 


sat. 


0.01831 


0.07124 


sat. 


0.01832 



Table 2. Characteristic times for two qubit independent decoher- 
ence for d = 3 dimensional density of states of field. i a ^ j a , 
ib + jb- 

We now analyze how the above results are af- 
fected when we consider the situation of collec- 
tive decoherence, as given by Eq. (16). The 
reduced density matrix for the two qubit system 
now reads: pf aJa ^ bJb (t) = exp {iS d {t) [( YL=a - 

Pfaja i b i h ^ n P ar ticular, we find: 



(i) ia = ja, ib ^ jb- 



Where /iia.idi f 

an d Joo.io = /li.oi 
ence rates are 



2i a (ib - jb)- Hence / 00 01 = / 11 10 = 1, 
-1. The corresponding decoher- 



pf =i (t) sa exp [ - Ti (i, T) ± id (uj c t + arctan(w c i))] p± =i (0) 



^= 3 W = ex P 



ic 3 2oj c t 



r 3 (*,r)± 

sin(2 arctan(o; c t)) 



1 + LOH 2 



PtM (26) 



for the Ohmic environment and the d = 3— dimensional 
density of states, respectively. 

(ii) i a = ja, ib = jb- 

Pf \ i i W = P? i i i as ex P ec t e d, the populations 
remain unaffected. 

(iii) ia ^ ja, ib ^ jb- 

pZ^S*) = ex p [- r -(*; T ) fc U,HA]^..v 4 (p)' where 

h 'i a j a ,hh = (ia + ib-ja - jb) 2 ■ Then, h ww = h' 01fil = 4, 
and h w 01 = /i 01 10 = 0. Obviously, the corresponding 
decoherence rates are: 



p; =1 w = p; =1 (o) , 

p+ i (t)«exp[-4r 1 (i,T)]p+ i (0) 
for the d = 1 dimensional field, and 

=Pia(0), 

exp [-4r 3 (t,T)]p+,(0) 



(27) 



(28) 



for the d = 3 dimensional field. Hence, regarding the 
input states, the case of collective decoherence shows two 
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very well defined situations: 

(a) A set of input states that shows no decoherence at all, 
despite the fact that the qubits are interacting with the 
environment. This is because under the specific situation 
of collective coupling there is a set of initial qubit states 
that does not entangle with the bosonic field and hence 
the states preserve their coherence. These states are the 
so-called "coherence-preserving" states and, for the case 
of a 2— QR, the corresponding reduced density matrix el- 
ements are (10|Tr B {p s (t)}|01), aadj01\Tr B {p s (t)}\10). 
As studied in more detail in Ref. |13| (where relaxation 
effects are also included), this result can be used as an 
encoding strategy, where an arbitrary L— QR can be de- 
coupled from its environment provided that every sin- 
gle qubit of the register can be encoded using 2 qubits: 
e.g. using the simple encoding (though not the most ef- 
ficient one) |0) i—* 1 01) , and |1) i— ► |10). It turns out that 
this procedure is ensured even if relaxation effects are in- 
cluded in the Hamiltonian (1) [Q. Hence in the specific 
situation of collective decoherence one can find, for arbi- 
trary L, a decoherence-free subspace (DFS) C L £ 7i® L 
(the Hilbert space H = Hqr(^Hb) that does not get en- 
tangled with the environment, and hence the QR should 
evolve without decoherence. Besides quantum error cor- 
rection codes, this is currently one of the most outstand- 
ing results in the battle against decoherence [|l3| , partic- 
ularly because of its relevance to maintaining a coherent 
qubit memory in quantum information processing. 

(b) The other two input states have the decay of coher- 
ences (ll|Tr B {p s (i)}|00) and (00\Tr B {p s (t)}\ll): these 
give a situation of 'superdecoherence' |g(J, where the 
qubits are collectively entangled and hence these matrix 
elements give the fastest decay for the coherences. This 
means that whilst there is a subspace which is robust 
against decoherence as discussed above, the remaining 
part of the Hilbert space gets strongly entangled with 
the environment. This superdecoherence situation is il- 
lustrated in Fig. 6 for the case of reservoirs with one and 
three dimensional-density of states. 

The above process of calculating explicit results for the 
decay of any coherence associated with the coupling of a 
L— QR to a bosonic reservoir, for both types of coupling 
(independent and collective), can be carried out for any 
L > 2 using the general formulas Eqs. (12) and (16). In 
so doing, we can obtain an estimate of typical decoher- 
ence times for a QR with an arbitrary number of qubits. 
We should point out that if no schemes for controlling the 
errors induced by the decoherence phenomenon are used 
|pr|-p7|, Tdec establishes an upper bound to the duration 
of any reliable quantum computing process. 



IV. DISCUSSIONS AND CONCLUDING 
REMARKS 

We have revisited a model of decoherence based upon 
the one previously studied by Leggett et al. [^2j in con- 
nection with the tunneling problem in the presence of 
dissipation, and used later by Unruh and Palma et al. 
p2| , p0t for describing the decoherence process of a quan- 
tum register. We have presented here a more complete 
description of this latter problem, which provides the fol- 
lowing new results. The decoherence rates of the den- 
sity matrix elements are correctly derived, leading to new 
quantitative results for both independent and collective 
decoherence situations. As discussed here, if no error 
correcting/preventing strategies are used, this has impli- 
cations for an estimation of decoherence time-scales for 
which the quantum memory of a QR can be maintained 
in any reliable computation (we note, however, there has 
been a recent proposal by Beige et al. regarding the use 
of dissipation to remain and manipulate states within a 
DFS fujl ). Our results agree with those reported in |||] 
for the case of collective decoherence but they are differ- 
ent to the ones reported there for the case of independent 
decoherence: in Ref. it is argued that independent 
decoherence, as opposed to collective decoherence, is in- 
sensitive to the qubit input states. Here we have shown 
instead that both cases are very sensitive to the input 
states and that both of them show collective decay. 

In the specific situation of independent decoherence we 
have found a non-trivial behaviour in the decay of the off- 
diagonal elements of the reduced density matrix. Here 
the coherences experience an effective revival before they 
either vanish or saturate to a residual value. This be- 
haviour depends on the temperature of the environment 
and depends strongly on the strength of the coupling: the 
coherence dynamics are different in the weak and strong 
coupling regime. Also, there are important qualitative 
differences between the Ohmic d — 1 and the super- 
Ohmic environment d = 3, which is ultimately linked to 
the spectral density of the bath. In the former case the 
coherence is always lost, whilst in the latter the coherence 
generally saturates to a residual value which is fixed by 
the temperature and strength of the coupling and only 
vanishes in the high-TE. By contrast, in the case of collec- 
tive decoherence we have identified QR input states that 
allow the system to evolve in a decoherence-free fashion, 
the so-called coherence-preserving states. We note that 
DFS's do not exist in the specific situation of indepen- 
dent qubit decoherence. We also note that our attention 
has been drawn to a dynamical-algebraic description that 
unifies the currently known quantum errors stabilization 
procedures §7j (see also Refs. @(a), §(|). Within the 
framework of a system formed by a collection of L uncor- 
rected clusters of subsystems where each cluster fulfills 
the requirements of collective decoherence, Zanardi has 
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shown that noiseless subsystems can be built p7[] ■ 

From the point of view of complexity analysis (and 
without including any strategy for the stabilizing of quan- 
tum information), we should ask how the results reported 
in the present paper affect those of Ref. p0[ ]. We must 
identify the coherences that are destroyed more rapidly: 
from Eqs. (12) and (16) it is easy to see that the coher- 
ences with the fastest decay are given by the matrix ele- 
ments p {0n l . and p, ,. These off-diagonal elements 
decay as expf— T,i(t;T)f(L)], with 
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APPENDIX A: SINGLE QUBIT DECOHERENCE 



f(L) = L + 2 ^ - 3m)(in ~ 3n) cos(uit s ) , (29) 



a=l,n=2 



for the limit of independent qubit decoherence, and 
as exp[— L 2 Td(t; T)] for the collective decoherence case. 
Hence, it is clear that for both cases, the longer the QR 
coherence length, the faster the coherence decay. Despite 
the fact that the results of Palma et al. |^(| are not the 
same as ours, it turns out that both sets of results lead 
to the same unwelcome exponential increase of the error 
rate. We note that the result of Eq. (29) is in general 
different to the one reported in ^0|. We also note that 
the coherence decay for the case of collective decoherence 
coincides with that of |2(J only for the fastest off-diagonal 
element decay: if we consider different density matrix el- 
ements, the results of Ref. |20) no longer coincide with 
ours. If the information reported in our work is used for 
an estimation of the actual decoherence time associated 
with any given off-diagonal density matrix element (co- 
herence), the results are in general quite different from 
the ones reported in 

We have shown how a bosonic environment destroys 
the coherences of an arbitrary quantum register. In do- 
ing so, we have identified DFS states that are invariant 
under the coupling to such an environment. This re- 
sult could be of crucial importance for improving the 
efficiency of quantum algorithms, for example. We be- 
lieve that the engineering of DFS's will become intrinsic 
to the designs of future quantum computation architec- 
tures. There was a recent experimental demonstration 
of decoherence- free quantum memories |lj|,[l9| . This has 
been achieved for one qubit, by encoding it into the DFS 
of a pair of trapped 9 Be + ions: in this way, Kielpinski 
et al. have demonstrated the immunity of a DFS of two 
atoms to collective dephasing [Q. Prior to this experi- 
ment with trapped ions, Kwiat et al. demonstrated the 
robustness of a DFS for two photons to collective noise 
p9[ . Robust quantum memories seem therefore to be 
well on their way, both theoretically and experimentally, 
to overcoming the main obstacle to quantum information 
processing - namely, decoherence. 

Acknowledgements. We are grateful to L. Viola for 
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The decoherence rates for a single qubit coupled to a 
reservoir with d = 1, and d = 3 density of states are (Eq. 
(19)) 



Ti(i,T) = p 



dujt 



•j/u 



1 — cos(wi) 



c ( »t!) ( —) (Al) 



(— 

V2cj t 



r 3 (t,T) = p J du we^ /Wc [l-cos(wt)]coth(^-). 

(A2) 

In Section III A we gave the analytic solutions to these 
integrals. However, we did not perform a full analysis 
of those results. We start here by recalling that the so- 
lution found for the integral (Al) was an approximate 
one, valid only for the low-TE (oj t << u) c ): a general 
solution to this integral requires numerical integration. 
The second integral was solved analytically for any tem- 
perature value T (without making any approximation). 
In these calculations note that the constant coupling ay 
changes its units with the dimensionality of the field: 
[aj = [(eV) 2 s 2 ], [a 3 ] =[(eV) 2 s 4 ], etc. 

We first analyze the case d — 1. In the low-TE, Eq. 
(20) leads to the identification of three main regimes for 
the decay of the coherences: (a) a "quiet" regime, for 
which t < t c , and Ti(t,T) w ciw 2 i 2 /2; (b) a "quantum" 
regime, where r c < t < r T , and Ti(t,T) w ciln(w c i); 
and (c) a "thermal" regime, for which t >> t t , and 
Ti(t,T) « 2ciOJ T t. These regimes have also been dis- 
cussed in Refs. [£0|,[32| and can be easily identified in Fig. 
7 for several different temperatures. In Fig. 7 (i) we have 
plotted Eq. (19) as a function of ui c t for several differ- 
ent temperatures and for d = 1. Since the decoherence 
effects arising from thermal noise can be separated from 
the ones due to quantum vacuum fluctuations, we have 
also plotted these partial contributions in order to see 
their effects over the time-scales involved in the decoher- 
ence of the single qubit (Eq. (Al)). It can be seen that 
for a given value of the temperature parameter 9, a char- 
acteristic time for which we start observing deviation of 
coherence from unity is determined by the shortest of the 
two time-scales r c and r T , and that this value is increased 
when decreasing the temperature T, as expected. 
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From Fig. 7 (i.c) it can clearly be seen that at low tem- 
peratures the quantum vacuum fluctuations play the ma- 
jor role in eroding the qubit coherence whilst the contri- 
bution due to thermal fluctuations plays a minor role. 
From this plot we can see the three main regimes indi- 
cated above: a quiet (t < r c ), a quantum (r c < t < t t ) 
and a thermal (t > t t ) regime. In this limit of low-TE, 
r c is the characteristic time that signals the departure of 
coherence from unity. Here the qubit dynamics shows a 
competition between contributions arising from vacuum 
and thermal fluctuations: even at thermal time scales, 
the contribution to the decoherence due to vacuum fluc- 
tuations remains important. 

In the case of the high-TE, the decay due to thermal 
noise (see dashed line in Fig. 7 (i.a)) becomes more im- 
portant than the vacuum fluctuation contribution f35| l 
and the start of the decoherence process is ruled by t t . 
Similar conclusions can be obtained from Fig. 7 (ii), 
where the decay of coherence has been plotted as a func- 
tion of time but in units of the thermal frequency oj t for 
several different temperatures. 

We can see from Table 3 and Fig. 7 that t decay is 
comparable to r c for the high-TE, and to r T for the low- 
TE. Indeed, if we assume cu c to be the Debye cut-off (cu c ~ 
10 13 s _1 ), we obtain from Table 3 that for c\ = 0.25, and 
T = 1 mK, the decoherence process starts at Tdec ~ 41.9 
fs and lasts for t decay ~ 27.4 ns (for d — 1). Here cu T — 
1.3 x 10 11 T » 1.3 x 10 8 s -1 , hence r T - 8 ns is of the 
same order of magnitude as tdecay For the super-Ohmic 
environment d = 3 (03 and T as before) we obtain Tdec ~ 
16.8 fs. In this case the coherences saturate to a residual 
value, as we discuss below (see Table 3 and Fig. 8 (ii)). 
From Table 3 we can also see the effects of the high-TE: 
the qubit decoheres several orders of magnitude faster 
than in the low-TE case. For example, for T — 100 K, 
ksT > hojc (hojc = 6.58 meV), and we obtain Tdec ~ 18.1 
fs, and tdecay ~ 0.6 ps (ci = 0.25, d = 1). A similar 
behaviour is also observed in the d = 3 case. 



d 




UJ T /iJc 


I^JcTdec 


e -ri(r dcc ) 


LO c tf 


e -r 4 (V) 






0.25 


10~ b 


0.418831 


0.98 


273950.34 


0.01 






0.25 


1.0 


0.181611 


0.98 


6.39891 


0.01 


d = 


1 


0.1 


10~ 5 


0.705612 


0.98 


1153307.91 


0.01 






0.1 


1.0 


0.291365 


0.98 


15.19703 


0.01 






0.01 


10~ 5 


7.47367 


0.98 


14346140.39 


0.01 






0.01 


1.0 


1.09604 


0.98 


147.12606 


0.01 






0.25 


10~ b 


0.167969 


0.98 


saturates 


0.778801 






0.25 


1.0 


0.154762 


0.98 


saturates 


0.564132 






0.25 


10 2 


0.020104 


0.98 


0.318417 


0.01 


d = 


3 


0.1 


KT 5 


0.275766 


0.98 


saturates 


0.904837 






0.1 


1.0 


0.251550 


0.98 


saturates 


0.795339 






0.1 


10 2 


0.031791 


0.98 


0.546769 


0.01 






0.01 


10 2 


0.101012 


0.98 


saturates 


0.135331 



In Table 3 we can appreciate the effect of the coupling 
strength on the decoherence time scales. Let c\ — 0.01 
(d = 1), hence for i) T = 1 mK, Td ec « 0.75 ps and 
tdecay ~ 1.4 for ii) T = 100 K, T dec ~ 0.11 ps and 
tdecay ~ 14.6 ps. Hence the weaker the coupling between 
the qubit and the environment, the longer the decoher- 
ence times and the slower the duration of the decoher- 
ence process. This result also holds for the case d = 3, 
as can be seen from Table 3. All of the above analy- 
sis concerning the different regimes for the decay of the 
coherences presented here is explicitly illustrated in the 
three-dimensional plots of Fig. 8. 

Next, we analyze the decoherence behaviour of the sin- 
gle qubit when coupling to the super-Ohmic d — 3 reser- 
voir. The corresponding decoherence function is given by 
Eq. (21). As can be seen from Fig. 8 (ii) and Table 3, this 
case shows an interesting behaviour for the coherence de- 
cay: once the end of the 'quiet' regime has been reached, 
the coherences decay to either zero, as in the case of Fig. 
8 (i) , or saturate to a particular value determined by the 
temperature parameter lu t /uj c . Here we can identify the 
particular temperature value for which no saturation oc- 
curs at all and the expected decoherence takes place. In 
Table 3 we give some saturation values for different tem- 
peratures and coupling strengths C3. For C3 = 0.1, and 
0.25 the temperature value for which any residual coher- 
ence vanishes falls in the interval 10 < 9 < 100 (high- 
TE). Apparently these residual values shown in Fig. 8 
(ii) and Table 3 vanish when additional frequency modes 
associated with the three-dimensionality of the field are 
taken into account (^0|. Even if this is not the case, 
surely the effects of relaxation mechanisms would mark 
the extent of these residual values. 



APPENDIX B: TIME EVOLUTION AND THE 
REDUCED DENSITY MATRIX 

In this Appendix we give details of the main steps fol- 
lowed in the calculation of the reduced density matrix 
given by Eq. (15). 



1. The time evolution operator Ui(t) 



Since Eq. (3) gives H t (t) = £ J?(tfe*V&t 



g*n e i " fc *5 fc ) ) we have that 



UM = Texp [- - j( ^(e^'j* b{ + 

e-^*'4 fc b h )dt'], (Bl) 



Table 3. Single qubit decoherence times for different temperatures, 
and coupling strength c; (t = 1,3); for d = 1 (Ohmic), and d = 3 
(super-Ohmic) dimensional density of states of the field. 



where we have introduced the shorthand notation = 
Y.n e ~ lk ' VnJ z- Hence, we can rewrite Eq. (Bl) as @: 
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U I (t)=e X p[j29 k ^ k (t)J k z bl] 

fc 

cxp (-J2g k ,^ k ,(t')Jz bl)g k ,P z k b k 



■ e -™ k t' x 



expfcwuAW? b l)} ■ (B2) 
fc' 

It is easy to show that the calculation of the product 
given by the last two lines Eq. (B2) gives the result 

g k Jt k [b h +g k ^ k (t)J*} . (B3) 

Hence, the following expression for U z (t) arises 

U r (t) = e^ h k k e ^ k k x 

exp { £ 5fc (t) J* &t - ^ (<) J]* 6 fe ] } 

fc 

cxp i \d* J z\ 



2 w 2 



(B4) 



where we have used the result e A+B = e A e B c ^ A ' B ^ 2 , 
which holds for any pair of operators A, B that satisfy 
[A, [A, B]]=0= [B, [A, B]] (as in the case of Eq. (B4)). 
It is straightforward to see that Eq. (B4) gives the final 
result 



_ , . r v ^ i 1 2 u. t — sinfo;. t) T+ . 



i jh 

Z J Z 



exp[^2{At(t)-A k (t)} 

fc 

where A k (t) = g* k <p* {t)J\ h b k . 

K k 



2. The reduced density matrix of a L — qubit register 

We start by using the result for C/j(i) in or- 
der to calculate the decay of the coherences, i.e. 
Tr B [ P B (0)U}^Ht)Uj-' l ^(t)], with Uj in} (t) as defined 
in Eq. (10). In so doing, we first compute the operator 
algebra for the product U^ n ^ (t)Uj 1 "^ (i) by taking into 
account the expression (B5). The result gives 



exp 



2 uj k t- sm(u) k t) 



k ■ r m „J exp |ff fc ^ fc (*)| 



^^(z m Z n jrnjn^ COS AC 

^ Wn sin fc • r m „ exp (cr fc &t _ cr*& fc ^ , (B6) 

m.n fc 

where we have set a k = g h <fu> k (*) X) m (*m - 3m)e~ lkrm . 
From the above equation note that the first two expo- 
nential terms commute, hence we only have to take the 
trace over the third term. By doing this (see e.g. Ref. 
p9|), we obtain the result 



ri exp 



9 k 



2 1 — cos(ui k t) , f TlUl,_ 



coth 



2k a T 



^ X ^m Jrn)(^n Jn) COS • 7*^ 

from where Eq. (11) arises immediately. 



(B7) 
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Figure Captions 

Figure 1. Two qubit "independent decoherence" due to 
the coupling to a reservoir of the Ohmic type (d = 1) as a 
function of time t and the transit time t 3 , for the input states 
associated with Ff(t,T) (i a / j a , it 7^ jb) ■ Here c\ = 0.25, 
and (i) 9 = 10~ 3 , (ii) 10°, and (iii) 10 2 . if (t,T), with i= 1,3 
are defined in the text. 

Figure 2. Two qubit "independent decoherence" caused 
by the coupling to an 'Ohmic environment' as a function of 
times t and t a , for the input states associated with T]~ (t, T) 
{ia + ja, i b + j b ). ci = 0.25, and (i) 9 = KT 3 , (ii) 10°, and 

(iii) 10 2 . 

Figure 3. Two qubit "independent decoherence" due to the 
super-Ohmic environment d — 3 (Eq. (25)) as a function of 
times t and t s , for the input states associated with Ff(t,T). 
c 3 = 0.25, and (i) 9 = 10 -3 , (ii) 10°, and (iii) 10 2 . 

Figure 4. Two qubit "independent decoherence" due to the 
super-Ohmic environment d — 3 as a function of times t and 
t s for the input states associated with T^(t,T). C3 = 0.25, 
and (i) 9 = 10" 3 , (ii) 10°, and (iii) 10 2 . 

Figure 5. Two qubit "independent decoherence". d = 3, 
coupling strength C3 = 0.01, and (i), (iii) 9 = 10 -3 , and (ii), 

(iv) 6» = 10 2 . 

Figure 6. Two qubit "collective decoherence" for (i) d = 
1 'Ohmic environment' (Eq. (27)), and (ii) d — 3 'super- 
Ohmic environment' (Eq. (28)), as a function of time and the 
temperature 9 = uj t /lu c . c\ = C3 = 0.25, and i a 7^ ja, and 
i b / jb- r+(t,T) (d = 1,3) is defined using Eqs. (27) and 
(28). 

Figure 7. (i) Decoherence of a single qubit for an 'Ohmic 
environment' as a function of t (in units of ui c ). The con- 
tributions arising from the separate integration of thermal 
(exp[— r T (t)]) and vacuum (exp[— T v (i)]) fluctuations are 
shown as dotted curves, ci = 0.25, (a) 9 = lo t /lo c = 1, 
(b) 10~ 2 , (c) 10" 5 . If lo c is the Debye cutoff, 9 « 10~ 2 T 
(see text): the decoherence shown corresponds to T = 100 K, 
T — 1 K, and T = 1 mK, respectively, (ii) Coherence decay 
for (a) 9 = 10" 5 , (b) 10" 2 , (c) 10 2 . ci = 0.25. Here time is 
given in units of the thermal frequency ui T = k B T/h. 

Figure 8. Decoherence of a single qubit for (i) d — 1, and 
(ii) d — 3, as a function of time (in units of the cut-off uo c ) 
and3 the temperature parameter 9 = uo T /uo c . ci = C3 = 0.25. 
If uj c is the Debye cut-off, the range of coherence decay goes 
from a few mK up to (a) 10 4 K (plot (i)) and (b) 1.5 x 10 3 K 
(plot (ii)). 



13 



This figure "figl.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/quant-ph/0105029v2 



This figure "fig2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/quant-ph/0105029v2 



This figure "fig3.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/quant-ph/0105029v2 



This figure "fig4.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/quant-ph/0105029v2 



This figure "fig5.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/quant-ph/0105029v2 



This figure "fig6.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/quant-ph/0105029v2 



This figure "fig7.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/quant-ph/0105029v2 



This figure "fig8.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/quant-ph/0105029v2 



